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PHASE ORDERING WITH AUTOMATIC WINDOW SELECTION fPAWS^ 
A NOVEL MOTION RESISTANT TECHNIQUE 
FOR 3D CORONARY IMAGING 

INTRODUCTION 

The development of navigator echo techniques has enabled the implementation of 
several algorithms to reduce motion effects from an image during acquisition (1-11). 
Early gating techniques, based on the acceptance/ rejection algorithm (2), were 
hindered byjhe loss in scan efficiency which resulted from the narrow acceptance 
windows required to effectively reduce artefacts and the variability of subjects' 
respiration. The use of phase encode ordering (8,9) and weighting (10,11) have been 
shown to improve image quality over conventional methods. The combination of 
phase ordering with a navigator acceptance window allowed larger acceptance 
windows to be used without compromising image quality whilst improving scan 
efficiency (9). Unfortunately, the effectiveness of all navigator acceptance window 
techniques is dependent on the choice of the acceptance window. Thus, during the 
scan, a change in the respiratory pattern may cause a window previously accepting 



data around the end-expiration position to be out of the new respiratory range. To 
overcome this limitation, more sophisticated gating techniques such as the 
Diminishing Variance Algorithm (DVA) have been suggested which adaptively 
select the most appropriate acceptance window depending on several factors such as 
the subjects' respiratory pattern, image data acquired so far, and scam time (4). The 
combination of DVA with a phase ordering or a weighting approach, although ideal, 
is not suited to changes in breathing patterns. The decision-taking is based upon a 
compromise between the repeated acquisition of profiles within an acceptable 
measuring~time and termination of the scan because the achieved image quality is 



already close to the desired one (11). Although these techniques attempt to be more 
adaptable to changes in respiratory patterns the use of an acceptance window is 
inherent as a decision to accept or reject data must still be made. This is therefore 
prohibitive if the breathing pattern changes significantly during the acquisition as 
respiratory positions which were previously discarded axe now desired and the 
previous rejection has resulted in a loss of scan efficiency. Although data can be 
acquired when the respiratory position is out of the required range, in case the 
diaphragm position changes during the scan, it is unpredictable how much will be 
useful when-the acceptance window is moved. The use of a weighted acceptance 
window further complicates this. 

The purpose of this work is to provide a technique which is resistant to 
changes in breathing whilst allowing the use of phase ordering to provide effective 
motion artefact reduction in the shortest time. Each position is regarded as equally 
relevant to the final image. The operator simply decides on the size of the 
acceptance window. 
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Viewed from one aspect the invention provides apparatus for magnetic 
resonance imaging a target object subject to periodic motion, said apparatus 
comprising: 

a first sensor for detecting when said target object is at a predetermined state 
within said periodic motion; 

a second sensor for detecting a signal indicative of a position of said target 
object relative to said apparatus; 

a magnetic resonance imaging scanner responsive to said sensor detecting said 
target object to be in said predetermined state for exciting said target object and 
recovering at least one line of imaging data in k-space; 

classifying logic for classifying said lines of imaging data into one of plurality 
of groups of lines of imaging data in dependence upon said position detected by said 
second sensor as said target object was excited, each group of lines corresponding to 
one of a plurality of contiguous ranges of position; and 

scan terminating logic for detecting when two or more groups of lines 
corresponding to contiguous ranges of position together contain a set of lines spanning 
k-space from which an image can be derived. 

It will be appreciated that the first sensor may be an ECG sensor and the second 
sensor may be an MR sensor detecting diaphram position. 




METHOD 

Theoretical Background 

The proposed technique, Phase-ordering with Automatic Window Selection (PAWS), 
uses a multi-level approach where no acceptance window is specified. Instead, the 
initial position of the diaphragm is taken as the reference and all further diaphragm 
positions are given an index position which is the displacement from this reference. 
Each index position is allotted a region of k-space, hereafter called a "bin", as shown 
in Figure 1(a). The k-space region is filled according to the direction of the arrows. 
Data acquisition for positions on the edges of k-space is straightforward, phase 
encode iine§~are acquired sequentially towards the centre of k-space. For positions 
which start in the centre of k-space, however, the data acquisition scheme is more 
complicated. Each diaphragm position can be used in one of a combination of 3 
separate bins. For example, for index number 4, the three possible combinations of 
bins which can complete k-space are: 2-3-4, 3-4-5 and 4-5-6. All three combinations 
are checked to determine which is closest to completion and this information is used 
to determine whether the right or left side of k-space is acquired. The side which has 
most unfilled k y lines remaining is acquired. This is demonstrated in Figure 1(b). 
For the case where the diaphragm is currendy at position 4, the window nearest to 
completion at this moment in time is 4-5-6. As there are more unfilled phase encode 
lines remaining between 4-6 than 4-5, the phase encode line toward diaphragm 
position 6 is acquired. Image acquisition is complete once the whole of k-space has 
been acquired by 3 contiguous bins. This is demonstrated in Figure 2. The final 
image is acquired in the shortest possible time and the phase encode lines have been 
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ordered thus minimising motion artefacts. In this example, the final range of motion 
is only 3mm as a bin size of 1 has been used with a navigator resolution of 1mm. 

Scan time can further be reduced in two ways: the bin size can be increased, 
thus allowing more than one index position to map to each individual bin, or the 
navigator resolution can be decreased so that a wider range of positions map to the 
same index position. However, image quality may be compromised as the phase 
ordering within each individual bin is less well defined. Figure 3 illustrates the final 
displacement graph for an image acquired with a bin size of 2 and navigator 
resolution of 1mm. Although the benefits of using phase ordering aire reduced with 
an increased^ bin size or lower navigator resolution, within each individual bin the 
displacement between any phase encode lines is limited to the bin size. Therefore, 
although a bin size of 1 and a navigator resolution of 2mm would produce a 
displacement graph as shown in Figure 2b, the range of motion in each bin is similar 
to that in Figure 3. 

To limit motion in the centre of k-space when the bin size is greater than one, 
a weighted bin scheme can be employed whereby an increased bin size is used for 
the two outer bins whilst the central bin size remains one. A typical displacement 
graph is shown in Figure 4. In this example, a bin size of 2 is used for the outer bins 
with a navigator resolution of 1mm for all the bins which gives a window size of 
5mm. 

In vitro Studies 

Studies were performed with a whole-body 0.5T magnet with a reduced 
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diameter (53cm LD.) gradient set giving 20mT/ m maximum gradient strength, and a 
slew rate of 60mT/m/ms. An MR imaging console (Surrey Medical Imaging 
Systems) with additional hardware and software was used to generate and drive the 
gradient and radio-frequency waveforms and to receive and reconstruct the image 
data. The ordering method was implemented for a fat-suppressed 3D velocity 
compensated gradient-echo pulse sequence (TE=6ms). 8 k z lines were acquired with 
a TR of 12ms, resulting in a total acquisition period of 96ms. Each k y was acquired 
after an interval of 800ms. 256 k y lines were acquired with a FOV of 300mm and a 
slab thickness of 2 0 mm resultin g i n a voxe l size of 1.2irunxl.2mmx2.5'mm. 

To compare the effects of the different ordering methods in vitro, complete 
data sets were acquired of a 4mm diameter tube. Multiple data sets were acquired 
with the tube moved by steps of 1mm between each. The total range of positions 
was 10mm. This allowed simulated data sets to be constructed given a series of 
motion measurements by selecting appropriate ky lines from each data set. The tube 
was filled with copper sulphate solution (0.5mM) and curved to represent a coronary 
artery. 

The PAWS technique was compared with the DVA. The DVA method allows 
a scan to be" terminated after a specified time or when the desired range of motion 
has been achieved. Therefore the methods can be compared to see how long each 
takes to acquire an image in a specified acceptance window and how effective 
respiratory artefact reduction is after a defined time. 

The techniques were compared with acceptance window sizes of 3mm 
(PAWS: bin size=l), 6mm (PAWS: bin size = 2) and 5mm (PAWS: bin size = 2 for 



7 



# # 

i =1 for central bin). A navigator resoluti 




outer bins, bin size =1 for central bin). A navigator resolution of 1mm was used in 
all cases. 

For scan efficiency, respiratory profiles acquired previously from 15 subjects 
were fed into a simulation program which calculated the scan time and diaphragm 
displacement graph through k-space for each technique, as well as the shortest seem 
time possible for each trace with the required acceptance window. 

For image quality, the DVA technique was used to construct an image in the 
same time as that taken by the PAWS technique. The displacement graphs through 
k-space produced by the simulation program were used to compile motion data sets. 
Profiles were taken through the images and, for image analysis, measurements were 
made of the peak signal intensity of the tube, the full width at half maximum 
(FWHM), the area surrounded by the pixels of the half maximum intensity (area at 
half maximum [AHM]) and the toted area (A). FWHM was used to measure blurring 
and (A- AHM) indicated the level of ghosting in the images. 
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RESULTS AND DISCUSSION 

Table 1 shows the number of cardiac cycles required to complete a scan for 
PAWS and the DVA method using the respiratory traces when acquiring an image 
of 128 PE lines. The shortest scan time possible is also displayed. The average 
percentage difference from the shortest scan time was only 0.5% (STD=0.7) for the 
PAWS approach in comparison to 15.5% (STD=26.3) for the DVA approach. The 
average percentage difference from the shortest scan time using both a weighted bin 
and a larger bin size are also shown in Table 2. Again, the PAWS approach achieves 
a shorter scan time with an average percentage difference of 4.6% and 4.7% 
respectively-Compared to 13.0% and 16.3% for the DVA approach. 

Figure 5 demonstrates the importance of navigator window selection and the 
implications for scan efficiency. A typical respiratory trace with the corresponding 
3mm acceptance window is shown. In this example, the PAWS technique (a) 
achieves optimal scan efficiency whilst the DVA technique (b) requires 47 extra 
cardiac cycles, a percentage difference of 18.1% from the optimal scan efficiency, to 
complete the scan. As explained later, the two methods have opted for slightly 
different acceptance windows. The DVA method includes the extreme end 
expiratory position in the acceptance window whilst the PAWS technique discards 
this position, instead opting to shift the acceptance window by 1mm to allow 
optimal scan efficiency. However, just a 1mm shift has caused a large reduction in 
scan efficiency. 

Figure 6 demonstrates the processes involved in selecting the final acceptance 
window. After collecting 128 cardiac cycles of data, the DVA method finds that the 
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most frequent diapru^n position is 287. The DVA tech^^ue will therefore attempt 
to narrow the range of motion around this position. The shortest scan time in this 
example is 260 cardiac cycles as this is the first point at which 128 diaphragm 
positions have been recorded within a 3mm window. As the PAWS technique has 
acquired everything, it detects that there are 3 contiguous positions which can be 
used to fill the whole of k-space and therefore terminates, using data in the range 
287-289. However, as the DVA technique is attempting to limit motion around the 
position 287, it is discarding data when the diaphragm is at position 289, even 
though we now know that the scan would finish with a greater scan' efficiency if the 
DVA technique had this position in its acceptance window. The scan must continue 
therefore, until 128 phase encoding lines have been acquired within the acceptance 
window 286-288. Although this is a suitable acceptance window, by studying the 
respiratory trace we can see in retrospect that better decisions could be made. 

This example highlights a drawback in the DVA technique. As it attempts to 
minimise motion around the most frequent position, it is assuming that this position 
will lie in the centre of the acceptance range. However, if this position is at end 
expiration, for example, which occurs often, then it is likely that it will occur towards 
one end of the acceptance window. Therefore, positions further away from end 
expiration, which should be included, are discarded as they are further from the 
maximum than other less frequent positions. This is more likely to occur with 
larger acceptance windows and could be a factor in the reduced efficiency of the 
DVA technique when the window size is increased. 
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The PAWS technique is also less efficient as the acceptance window is 
increased and the principal reason for this lies in the grouping of diaphragm 
positions when the bin size is greater than one. The initial diaphragm position 
which is recorded is set to index 0 and each successive diaphragm position is given 
an index relative to this. In cases where the bin size is greater than one, diaphragm 
positions are grouped together in what may eventually not be the most optimal 
way. This is demonstrated in Figure 7. For a 5mm window, the most efficient scan 
time for the acquisition of 128 phase encode lines would have occurred after 190 
cardiac cycles in the window 287-291. However, because of the 'way PAWS has 
grouped the diaphragm positions, 286 and 287 are paired together and the window 
287-291 is not possible. The acquisition must therefore continue and completes after 
217 cardiac cycles. 

If a window combination is therefore unavailable to PAWS, there is potential 
for loss of scan efficiency given the scan parameters. However, as we mentioned 
earlier, there are two possible ways of achieving the same window sizes: by altering 
the bin size and by altering the navigator resolution. Therefore, if we are using a bin 
size of 2 and navigator resolution of 1mm, PAWS may not perform in the shortest 
time possible as the most efficient window is not available. If these parameters are 
changed, however, and a bin size of 1 is used with a lower navigator resolution of 
2mm, the window size remains the same but PAWS can now complete in the 
shortest time possible for the scan given the parameters. The use of a reduced 
navigator resolution for the 6mm window reduced the average difference from the 
shortest scan time from 4.7% to 2.6%. It is important to note however that although 
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an improvement wastound by using the reduced navigator resolution as opposed to 
a larger bin size / this method is also grouping diaphragm positions, the difference 
being that the PAWS algorithm can now complete with the same window as that 
chosen by the simulation program when calculating the shortest scan time possible. 
Depending on the grouping, the shortest scan time possible may change from one set 
of parameters to the next. However, given the infinite variations of possible 
parameters, as the PAWS technique is more efficient given a single bin size, we feel 
that the lowest navigator resolution possible should be used. This technique also 
aided the DVA approach. In one particular case (subject 1) the PVA method 
selected a window around inspiration as the most frequent position occurred here. 
However, grouping navigator positions by way of reducing the resolution of the 
navigator allowed DVA to focus on a larger more frequent region and allowed the 
scan to terminate in a much shorter scan time with a saving of 203 cardiac cycles. 

Figure 8 demonstrates the problems encountered when attempting to acquire 
data in a situation where the breathing pattern is extremely variable. The 
accumulated diaphragm position histograms (a) are shown alongside the differential 
histograms (b) and the corresponding breathing pattern is shown in Figure 9. This is 
an actual breathing pattern and demonstrates the inconsistent nature of breathing 
often observed during an acquisition. In this example, 256 views are acquired to 
demonstrate the changes in breathing and implications for scan efficiency and image 
quality during a long acquisition. As we can see, after 256 cardiac cycles, diaphragm 
position 276 is the most frequent and, therefore, the DVA technique attempts to 
minimise motion around this position. However, during the next 100 cardiac cycles, 



this position only occurs twice as the breathing pattern has changed significantly. If 
we continue to monitor the differential histogram we see that the breathing pattern 
differs greatly from one period of 100 cardiac cycles to the next and eventually 
reverts back to a pattern more similar to the original. Therefore, there would be little 
advantage in altering the acceptance window for the DVA technique as we would 
simply encounter the problem of a change in the breathing during each successive 
period. With such a varying breathing pattern, a technique which bases it decisions 
on the breathing pattern so far to predict what will happen next is unlikely to be 
particularly successful and it is in such cases that the PAWS technique is most useful. 
As PAWS is_ ordering and accepting all the data, it simply completes once the 
appropriate number of phase encoding lines have been acquired within contiguous 
bins. Therefore, regardless of the breathing pattern, the PAWS algorithm will 
complete in, or close to, the shortest possible time. In this example, PAWS completes 
in 819 cardiac cycles whilst DVA requires 1025 cardiac cycles to complete. 

Figure 9 also illustrates a problem of all techniques which use the frequency 
information to select the most suitable window. Both DVA and PAWS select a 
window which is towards inspiration (DVA: 274-278; PAWS: 270-274). However, 
acquisitions^ made when the diaphragm is towards inspiration are more likely to 
suffer motion artefact as the diaphragm is moving with a greater velocity (12). 
Therefore, it is not only important to place the acceptance window around a 
frequently occurring diaphragm position, but this position should also preferably be 
around end-expiration. Again, the PAWS technique is ideally suited to this 
requirement as data from all respiratory positions has been acquired. Therefore, 
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encode lines have been comple 




3 contiguous bins, data 



acquisition can continue with the aim of completing all the phase encoding lines in 



an end-expiration window. As we will know how close to completion each window 
is, we can judge how long the additional scanning is likely to require and if data 
acquisition takes too long, then the images from the original window can be 
accepted by the user. In the example in Figures 8-9, a 5mm window could have been 
completed by PAWS around end-expiration after 919 cycles, more than 100 cardiac 
cycles before the DVA technique completed. 

It is also important to note that the window, [c], which has been chosen to be 
the most desirable in terms of image quality in Figure 9 is unlikely to have been 
selected were a decision to be formed after either an initial period of monitoring the 
diaphragm position or after analysis of scan efficiency as many of the positions occur 
much later in the scan and many do not occur even in the initial 256 cardiac cycles. 
Instead, it is- only in retrospect, having weighed up both scan efficiency and image 
quality costs, that this window is selected. This is the principal benefit of the PAWS 
technique. It is possible during data acquisition to analyse the subject's breathing 
pattern and the progress of acquisition for any window. As the order of data 
acquisition is the same, regardless of the final window we may choose, it is possible 
to make an infinite number of alterations to the final window with no loss in scan 
efficiency. This is extremely useful in situations where a window placed around 
end-expiration is now acquiring data in inspiration. 

As well as the desire to acquire images from around end-expiration, there is 
also a preference as to the final ordering of the phase encoding lines. In our previous 
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studies into the effects of different ordering patterns on coronary artery imaging, we 
found that the use of an asymmetric ordering pattern produced considerably better 
images than the use of a symmetric ordering pattern (13). Therefore, if the initial 
data set which completes produces a symmetric ordering pattern, we could balance 
up the expected gains in image quality and the extra time required to complete a 
more suitable window resulting in an asymmetric pattern. It is important to note, 
however, that both ordering patterns produced better images than those acquired 
without the use of ordering. 

Eigure_10_demonstrates-the-advantages-of--txsing--phase ordering-with the 

PAWS technique. In each case the phase ordered images have reduced ghosting and 
blurring. However, as we would expect, as the window size is increased the 

— artefacts-increase for both methods (Fig. lOb-d). The mean and standard deviation 
for ghosting and blurring measurements for each method are shown in Table 3 for 
all 15 subjects. The PAWS method was more effective at reducing both blurring and 
ghosting artefacts if the motion in the centre of k-space was limited (Fig. lOa-b). As 
we increase the bin size, the advantages of using ordering are diminished (Fig. 10c- 
d). However, the use of a weighted technique, whereby the motion in the centre of 
k-space is restricted whilst the outer regions are acquired within larger acceptance 
ranges, allows an effective compromise between image quality and scan efficiency 
(Fig. 10b). 

As the PAWS and DVA techniques select a window primarily on the basis of 
scan efficiency, it is possible that over successive acquisitions the position of the 
window could change. This could lead to misregistration problems between images. 
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However, it could be possible to combine the techniques with slab following to 
ensure a more consistent window for all images. Therefore, after the initial 
acquisition, the centre of the window is used as a reference position. We could use 
slab following to ensure that all further acquisitions occur at the same point. The use 
of phase ordering could also help overcome the small phase errors and pseudo- 
motion ghosting from stationary tissues which result from the use of slab following 
by ensuring that the relative errors between successive phase encode lines is 
minimal. 

Once the appropriate window has been selected from the multiple windows 
which have _been acquired, several partially complete data sets remain. It may 
therefore be possible to use a combination of data averaging and post-processing to 
utilise more of the data to further improve image quality. Also, more than one data 
set could be acquired in a single acquisition. 

Conclusion 

A method is introduced which is resistant to changes in breathing and allows images 
to be acquired in the shortest possible scan time with no requirement for operator 
interaction. The ability to integrate the use of phase ordering with a flexible window 
selection strategy further allows improvements in image quality. This method has 
proved to be effective both in improving scan efficiency and reducing respiratory 
motion artefacts in our in vitro studies. We believe this technique will aid in 
overcoming many of the current problems faced with navigator acceptance 
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techniques, providing an ordering technique robust against changes in breathing 
which has previously not been possible. 
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TABLES 





Least Cardiac 
Cycles Required 


DVA Cardiac 
Cycles 


PAWS Cardiac 
Cycles 


1 


559 


1103 


565 


2 


507 


671 


507 


3 


765 


765 


765 


4 


260 


307 


260 


5 


300 


300 


300 


6- 


-238 


268 


239 


7 


290 


290 


291 


8 


273 


273 


273 


9 


227 


227 


230 


10 


515 


520 


515 


11 


467 


467 


467 


12 


451 


497 


453 


13 


358 


420 


366 


14 


541 _ " 


774 


547 


15 


364 


364 


364 



Table 1, Number of cardiac cycles required to acquire 128 phase encode lines within 
a 3mm acceptance window. 
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* # 



Acceptance Window 


DVA 


PAWS 


3 (Bin 1-1-1) 


15.5% (STD 26.3) 


0.5% (STD 0.7) 


5 (Bin 2-1-2) 


13.0% (STD 17.9) 


4.6% (STD 4.8) 


6 (Bin 2-2-2) 


16.3% (STD 19.1) 


4.7% (STD 4.6) 



Table 2. Percentage difference between cardiac cycles required to acquire 128 phase 
encode lines using the DVA and PAWS technique compared to the optimum scan 
time. 
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Acceptance Window 


DVA 


PAWS 


3 (Bin 1-1-1) 


4.4% (STD 0.4) 


4.2% (STD 0.5) 


5 (Bin 2-1-2) 


4.7% (STD 0.5) 


5.3% (STD 1.2) 


6 (Bin 2-2-2) 


5.1% (STD 0.7) 


5.0% (STD 0.8) 



Table 3a, Mean blurring measurements for tube using traces acquired from 
subjects using PAWS and DVA. 
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# 



Acceptance Window 


DVA 


PAWS 


3 (Bin 1-1-1) 


434.7% (STD 44.8) 


330.9% (STD 12.8) 


5 (Bin 2-1-2) 


423.3% (STD 61.5) 


371.6% (STD 28.8) 


6 (Bin 2-2-2) 


458.6% (STD 67.2) 


438.9% (STD 65.8) 



Table 3b. Mean Ghosting measurements for tube using traces acquired from 15 
subjects using PAWS and DVA. 
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Legends to Figures 

1. Each index position is allotted a starting position in k-space, denoted by the solid 
circle. K-space is filled according to the direction of the arrows. The data 
acquisition scheme in the centre of k-space is more complex. In this example, for 
diaphragm position 4, the next phase encode line to the right of k-space would be 
chosen as window 4-5-6 is the closest to completion and there are more phase 
encode lines left to acquire between diaphragm positions 4 and 6. 

2. Data acquisition is terminated once the whole of k-space has been acquired by 3 
contiguous bins (a). The final displacement through k-space is ordered to reduce 
respiratory artefacts (b). 

3. If the bin size is increased, each individual bin is not ordered which could result 
in a loss of image quality. 

4. To limit motion in the centre of k-space, a weighted bins method can be 
employed. In this example, a bin size of 2 used for the two outer bins whilst the 
central bin size was one. 

5. The PAWS technique is more efficient than the DVA technique in this typical 
example. Notice the different windows each method uses. The DVA technique 
has included the extreme end expiration position whilst the PAWS technique has 
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discarded this position and shifted the acceptance window by 1mm to allow a 
more efficient scan time. 

6. After 128 cardiac cycles, the most frequent diaphragm position is 287. The DVA 
technique aims to narrow the range of motion around this position. The 
optimum scan time in this example is 260 cardiac cycles as it is only at this point 
that 128 diaphragm positions have been recorded within a 3mm window. The 
PAWS technique therefore terminates and the window is set to the range 287-289. 
The DVA technique must continue as it is attempting to limit motion around 287 
and as 289 is further away than 286 it is discarded. 

7. This example illustrates how the grouping of diaphragm positions may result in 
scan inefficiencies. A weighted scheme is used of 2 for the outer bins and 1 for 
the central bin to produce a 5mm acceptance window. Positions which have been 
grouped together are shown linked with a short bar in the graph above. In this 
example, 128 phase encode lines could be most quickly acquired within the 
window 287-291. However, as PAWS has grouped together positions 286-287 
and 289-290, data acquisition continues until enough lines have been acquired in 
the window 286-290. 

8. When the breathing pattern is very variable it is difficult to select an appropriate 
window. After acquiring 256 views, the DVA method attempts to minimise 
motion around 276 (a). However, for the next 100 cardiac cycles, the breathing 

i 
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pattern changes so that very few positions are registered around 276 (b). It 
would not be beneficial to change the acceptance window at this stage as the 
breathing pattern continues to change over successive periods and again shifts 
back to a state similar to its original. In this example, the DV A technique requires 
1025 cardiac cycles to complete whilst the PAWS technique completes in 819 
cardiac cycles. The 'ideal' end-expiration window, which has been 
retrospectively selected, completes after 919 cardiac cycles. 

9. This extremely variable breathing pattern makes selection of an acceptance 
window- very difficult. The acceptance window selected by the DVA technique 
[a] and PAWS [b] is highlighted along with a desirable window around end- 
expiration [c]. 

10. Phase ordering, alongside the improved scan efficiency of the PAWS technique, 
produces better images than the DVA technique in the same scan time. 
However, as the bin size is increased, particularly in the central region of k-space 
(c-d), the benefits of using phase ordering are reduced. 
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ABSTRACT 

Navigator acceptance imaging methods are hindered by the loss in scan efficiency 
which results from the changes in the breathing pattern of a subject over time. The 
Diminishing Variance Algorithm (DVA), which does not use a predefined 
acceptance window, is less influenced by such changes. The use of phase ordering 
and weighting techniques have been shown to significandy improve image quality 
over non-ordered window methods. However, the use of an acceptance window is 
inherent in all these techniques as a decision to accept or reject data must still be 
made. A technique is presented which is resistant to changes in breathing whilst 
allowing the'use of phase ordering to provide effective motion artefact reduction in 
an optimal time. The basic principle is described and illustrated for this automatic 
window selection technique and in vitro results are presented to demonstrate the 
feasibility of this method. 

• Coronary artery imaging 

• Respiratory Gating 

• Motion artefacts 

• Phase encode ordering 

[Figure 2] 




0 ky line +54 



Figure 1. Each index position is allotted a starting position in k- 
space, denoted by the solid circle. K-space is filled according to the 
direction of the arrows. The data acquisition scheme in the centre of 
k-space is more complex. In this example, for diaphragm position 4, 
the next phase encode line to the right of k-space would be chosen 
as window 4-5-6 is the closest to completion and there are more 
phase encode lines left to acquire between diaphragm positions 4 
and 6. " 




Figure 3. If the bin size is increased, each individual bin is not 
ordered which could result in a loss of image quality. 
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Figure 4. To limit motion in the centre of k-space, a weighted bins 
method can be employed. In this example, a bin size of 2 used for 
the two outer bins whilst the central bin size was one. 
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Figure 6, After 128 cardiac cycles, the most frequent diaphragm position is 287. The 
DVA technique aims to narrow the range of motion around this position. The 
optimum scan time in this example is 260 cardiac cycles as it is only at this point that 
128 diaphragm positions have been recorded within a 3mm window. The PAWS 
technique therefore terminates and the window is set to the range 287-289. The DVA 
technique must continue as it is attempting to limit motion.around 287 and as 289 is 
further away than 286 it is discarded. 
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Figure 7. This example illustrates how the grouping of diaphragm positions may 
result in scan inef ficienci es. A weighted scheme is used of 2 for the outer bins and 1 
for the central bin to produce a 5mm acceptance window. Positions which have been 
grouped together are shown linked with a short bar in the graph above. In this 
example, 128 phase encode lines could be most quickly acquired within the window 
287-291. However, as PAWS has grouped together positions 286-287 and 289-290, 
data acquisition continues until enough lines have been acquired in the window 286- 
290. 
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Figure 8. When the breathing pattern is very variable it is difficult to select an 
appropriate window. After acquiring 256 views, the DVA method attempts to minimise 
motion around 276 (a). However, for the next 100 cardiac cycles, the breathing pattern 
changes so that very few positions are registered around 276 (b). It would not be 
benef icial to change the acceptance window at this stage as the breathing pattern 
continues to change over successive periods and again shifts back to a state similar to its 
original. In this example, the DVA technique requires 1025 cardiac cycles to complete 
whilst the PAWS technique completes in 819 cardiac cycles. The 'ideal' end-expiration 
window, which has been retrospectively selected, completes after 919 cardiac cycles. 




Figure 9. This extremely variable breathing pattern makes selection of an acceptance 
window very difficult. The acceptance window selected by the DVA technique [a] 
and PAWS [b] is highlighted along with a desirable window around end-expiration 
[c]. 
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DVA 
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Blurring: 4.52 
Ghosting: 361 



(a) 



Blurring: 4.75 
Ghosting: 476 



(b) 



Blurring: 5.2 
Ghosting: 490 



(c) 



Blurring: 6.4 
Ghosting: 535 




3mm 



Blurring: 4.12 
Ghosting: 320 



5mm 



-Blurring: 3.92 
Ghosting: 355 



6mm 



Blurring: 5.1 
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Figure 10- Phase ordering, alongside the improved scan efficiency of the PAWS 
technique, produces better images than the DVA. technique in the same scan time. 
However, as the bin size is increased, particularly in the central region of k-space (c- 
d), the benefits of using phase ordering are reduced. 
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